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ABSTRACT 

It is generally accepted that Anomalous X-ray Pulsars (AXPs) and Soft Gamma-ray Repeaters 
(SGRs) are magnetars, i.e. neutron stars with extremely high surface magnetic fields {B > 
lO^"* G). The origin of these high magnetic fields is uncertain, but a popular hypothesis is that 
magnetars are born with an initial spin period not much exceeding the convective overturn 
time (~ 3 ms), which results in a powerful dynamo action, amplifying the seed magnetic field 
to > 10^^ G. Part of this rotation energy is then expected to power the supernova through 
rapid magnetic braking. It is therefore possible that magnetars creation is accompanied by 
supernovae that are an order of magnitude more energetic than normal supernovae, provided 
their initial spin period is ^^ 1 ms. However, we list here evidence that the explosion energies 
of these supernova remnants associated with AXPs and SGRs - Kes 73 (AXP IE 1841-045), 
CTB 109 (AXP 1E2259H-586) and N49 (SGR 0526-66) - are close to the canonical supernova 
explosion energy of 10^^ erg, suggesting an initial spin period of Po ?i 5 ms. 

We therefore do not find evidence that magnetars are formed from rapidly rotating proto- 
neutron stars, allowing for the possibility that they descend from stellar progenitor with high 
magnetic field cores, and we discuss the merits of both formation scenarios. 

In an appendix we describe the analysis of XMM-Newton observations of Kes 73 and 
N49 used to derive the explosion energies for these remnants. 

Key words: stars:neutron - stars: magnetic field - ISM: supernova remnants - 
ISM:individual:N49 - ISM:individual:Kes 73 - ISM:individual:CTB 109 



1 INTRODUCTION 



The notion that neutron stars exist with surface magnetic fields 
as high as 10^* — 10^^ G has become generally accepted over 
■ the last decade. The most spectacular manifestations of these 
' so-called magnetars are arguably the Soft Gamma-ray Repeaters 
(SGRs), which have received ample attention after the giant flare 
of SGR 1806-20 on December 27 2004 which had an energy 
of ~ 10*'' erg (e.g. iHurlev et al] l2005l) . Related to the class 
of SGRs are the Anomalous X-ray Pulsars (AXPs), which are 
less prone to flare, a lthough occasional flares have been observed 
JGavriil et alj|2002t iKaspi et a l."2003'). Like SGRs they are ra- 
dio quiet X-ray pulsars with spin periods in the range of 5-12 s. 
The timing properties of both SGRs and AXPs suggest that they 
are isolated neutron stars, whereas their spin down rate suggest 
that they have ma gnetar-like dipole magneti c fields of Bdip ~ 
10" " 10^^ G (see Woods & T hompson 2004, for a review). Other 
arguments in favour of the magnetar hypothesis can be found in 
iThompson & DuncarJ il993) . 

The apparent dichotomy between the radio quiet magnetars, 
on the one hand, and ordinary young radio pulsars, with B ~ 
lO^'^ — lO^'' G, on the other hand, is likely to originate from distinct 
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properties of the progenitor stars. One idea is that the high magnetic 
field of magnetars simply reflects the high magnetic field of their 
progenitor stars. Magnetic flux conservation iWoltieill964) implies 
that magnetars must then be the stellar remnants of stars with in- 
ternal magnetic fields of B > 10"* G, whereas radio pulsars must 
be the end products of stars with B ~ 10'^ G. Based on a st udy 
of magnetic white dwarfs, Ferrario & Wickramasinghg i200a) ar- 
gue that there is a wide spread in white dwarf progenitor magnetic 
fields, which, when extrapolated to the more massive progenitors 
of neutron stars can explain the existence of magnetars. 

A currently more popular idea for the origin of magnetars 
is that they are formed from proto-neutron stars with periods in 
the range Pi ~ 0.6 ms (the break-up limit) to 3 ms. Convec- 
tion and differential rotation in such rapidly rotating stars would 
give rise to an efficient a — Sl-dynamo, resu lting in magnetic field 
amplification on a time scale of < 10 s CDuncan & Thompsod 
1992; Thompson & Duncan 1993; Duncan & Thompson 1996J). 
In contrast, radio pulsars presumably start their life with initial 
spin periods of > 10 ms o r longer. As initially pointed out by 
Dun can & ThompsonI il992h . the idea that magnetars form from 
proton-neutron stars with P ~ 1 ms implies that supernova explo- 
sions that create magnetars are an order of magnitude more ener- 
getic than ordinary core collapse supernovae. The reason is that a 
neutron star spinning with P = 1 ms (fi = 6.3 x 10'' s~^) has a 
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Table 1. The explosion energies and ages of tlie supernova remnants from X-ray spectral analysis. 



SNR/Pulsar 



Distance 
kpc 



radius 
pc 



E 
(10^1 erg) 



?iH 



Mass 



SNR Age 
lO^yr 



Pulsar Age 
10=^ yr 



References 



Kes73/lE1841-()45 


7.0 


4.3 


0.5 ±0.3 


2.9 ±0.4 


29 ±4 


1.3 ±0.2 


4.3 


CTB109/1E2259+586 


3.0 


10 


0.7 ±0.3 


0.16 ±0.02 


97 ±23 


8.8 ±0.9 


220 


N49/SGR 0526-66 


50 


9.3 


1.3 ±0.3 


2.8 ±0.1 


320 ± 50 


6.3±1.0 


1.9 



This work 

Sasaki et al. (2004) 

This work 



Note that lSasaki et alji2004h find an energy of fl.g ± 0.7) X 10"^^ erg for CTB109, if they assume incomplete te mperature equilibration at th e 
shock front. Derived energies scale with distance d as d?'^ . Distances and pulsar ages (r = -^P/P) are taken from lWoods & Thompson! 120041) . 



rotational energy of Erot = ^I^'^ ^ 3 x W ^^jj^) erg, for 
a moment of inertia of / ~ 1.4 x 10*^ gcm^ ( Latti mer & PrakashI 
12001) . Due to the high magnetic field, the magnetar rapidly spins 
down during the explosion, thereby powering the supernova. An 
upper limit for the time scale in which most of the rotational energy 
is lost by a neutron star with B — lO^'' G is obtained by assuming 
magn etic braking in vacuum, r ~ 160 s. More detailed model- 
ing bv lThompson et alj J2004) indicates that the time scale for spin 
down of the proto-neutron star may be as short as 10-100 s. 

Given the implication of the millisecond proto-neutron star 
hypothesis that magnetars should be surrounded by hyperiiowa rem- 
nants, it is surprising that this idea has not directly been tested on 
supern ova remnants (SNRs) as sociated with AXPs and SGRs, al- 
though |Aiie^^Ho^3 IZQol) did calculate the effect of the hy- 
pernova/magnetar connection on the evolution of the SNR, and ap- 
plied this to known magnetar/SNR connections. Here we review 
what is known about the explosion energies of the supernova rem- 
nants Kes 73, N49 and CTB 109 associated with respectively the 
magnetar candidates AX? 1E1841-045, SGR 0526-66, and AXP 
1E2259-I-586. For CTB 109 and N49 explosion energy estimates 
have been given in the literature, but without discussing the possi- 
ble implications for the magnetar formation scenario. As far as we 
know, no explosion energy estimates exists for Kes 73, we therefore 
present our own energy estimate based on archival XMM-Newton 
observations. In addition we present an analysis of the explosion 
energy of N49, also based on unpublished XMM-Newton data. 



THE EXPLOSION ENERGIES OF SUPERNOVA 
REMNANTS ASSOCIATED WITH MAGNETARS 



medium JTruelove & McKeall99St) . The shock radius, rs, and ve- 
locity, Vs , of such an explosion evolves as 



iGaenslen l2004h reviewed the association of AXPs and SGRs with 
SNRs, taking into account the probability of random spatial coin- 
cidence given the extent of the SNR and the location of the AXP or 
SGR within the SNR. He concludes that magnetar candidates likely 
to be associated with SNRs (i.e. have a chance alignment probabil- 
ity < 1%) are: AXP 1E1841-045 with Kes 73 (G27.4-I-0.0), AXP 
1E2259-I-586 with CTB 109 (G109.1-1.0), SGR 0526-66 withN49, 
and the candidate AXP AX J 1845-045 with G29.6-^0.1. 

G29.6-I-0. 1, a shell-ty pe SNR, was detected in archival radio 
data iGaensler et al. 1999) following the discovery of AX J1845- 
045 I Gotthelf & Vasisht 1998). Unfortunately, the X-ray emission 
from the SNR is too weak ( Vasisht et al. 2000) and the distance too 
poorly constrained iGaensler et al.l9 99) to derive any meaningful 
energy estimate from the X-ray spectrum of the SNR. This leaves 
us with Kes 73, N49 and CTB 109 to test whether SNRs with mag- 
netars are an order of magnitude more energetic than regular SNRs. 

One can determine the energy of SNRs using the well known 
Sedov evolution of an adiabatic point explosion in a uniform 
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with po the interstellar medium density, E the explosion energy, 
and t the age of the SNR. The shock velocity can be determined 
from X-ray spectra, because for a monatomic gas and high Mach 
number shocks the following relation holds between post-shock 
plasma temperature and shock velocity: 



3 2 

kT = -—0.63mpVs, 
16 



(3) 



with rrip the proton mass. The factor 0.63 gives the average particle 
mass in units of the proton mass. It is smaller than the proton mass, 
because it includes electrons. So from measuring the temperature 
one can estimate Vs, which, combined with the measured radius, 
gives the age of the SNR (Eq. Q&|2}. The interstellar medium 
density can be estimated from the X-ray emission measure (EM 
= J UenndV). Having obtained a density estimate, the explosion 
energy can be calculated using Eq.Q 

We list the results of our analysis of XMM-Newton data in 
Table [H together w ith the parameters for CTB 109 reported by 
Sasa ki et alj J2004|). It is clear that these SNRs with magnetars do 
not have an order of magnitude higher explosion energies than the 
canonical supernova explosion energy of 10^^ erg, suggesting an 
initial period P > 5 ms. More details about our analysis data can 
be found in the appendix, which also includes a short discussion 
on the effects of non-equilibration of temperatures on the derived 
parameters, and on the applicability of the Sedov model for Kes 73. 
FigureQshows the XMM-Newton spectra of Kes 73 and N49 with 
best fit spectral models. 

For all three SNRs the distances are well determined. N49 
is located in the Large Magellanic Cloud, Kes 73 has a reliable 
distance estimate based on HI absorption measurements (d = 
6 - 7.5 kpc, Sanbonmatsu & Helfand 1992), and CTB 109 can- 
not be more distant than Cas A, which is situated at 2° from CTB 
109, has a higher absorption column, and is at a distance of 3.4 kpc 
JReed et alll995l) . Both TB109 and Cas A are likely to be associ- 
ated with the Perseus arm at ~ 3 kpc. 

As an aside, Cas A contain s a point source that is suspected 
to be a magnetar (see Pavlov et al. 2004, for a discussion), a sus- 
picion that has gained more credibility with the recent discovery 
of an infra-red lig ht echo from a putative SGR-like giant flare 
iKrause et all200 5). However, Cas A appears to have an explosion 
energy of ~ 2 x 10''^ erg, which is quite accurately known, as the 
expansion has been directly measured ( Vink 2004|, for a review). 

We like to conclude this section by pointing out that, although 



© 2005 RAS, MNRAS 000,[T]-?? 



Supernova remnant energetics and magnetars 3 





MB 


9IH 


1 


Sli'jil 




B 


1 






I^H^IHI 



O0.5 




Energy (keV) 



Energy (keV) 



Figure 1. XMM-Newton MOS12 spectra of Kes 73 (left) and N 49 (right). The sohd hnes represent the best fit Sedov models. See appendix for details. Chandra 
ACIS images (taken from the archive) are shown as insets. The RGB color coding for Kes 73 con'esponds to the energy ranges: 0.6-1.5 keV, 1.5-2.8 keV, 2.8-5 
keV. The image measures 5.1' X 5.1'. The N49 images consists of a mosaic, and the RGB colors correspond to 0.5-0.75 keV, 0.75-1.2 keV and 1.2-3 keV. The 
image measures 1.9' X 1.9', but due to the small window mode of the observations the southwestern comer of N49 is missing. AXP 1E1841-045 is the bright 
unresolved source in the center of Kes 73; SGR 0526-66 is the unresolved source near the northern edge of N49. For both images a linear brightness scaling 
was used, but saturating the unresolved sources in order to bring out the SNR emission. 



there are some uncertainties associated with determining the ex- 
plosion energies using a Sedov analysis (gradients in the medium 
in which the SNR evolves, conversion of part of the energy to cos- 
mic rays), these uncertainties apply equally to all other SNRs for 

which energies have been determined. Nevertheless, other SNRs 

I I I J , 

appear to have similar explosion energies (Huehes et al. 199b). In 

other words, the most remarkable feature of SNRs with magnetars 
is that they are unremarkable, i.e. they are indistinguishable from 
normal SNRs. 



3 DISCUSSION 

This study shows that in general magnetar formation is not ac- 
companied by a hypemova explosion. This does not necessarily 
preclude the idea that observed extra-galactic hypernovae, such 
as SN1998bw, are not by the formation of a rapidly spinning 
magnetar followed by rapid magnetic braking iNakamura 1998; 
iThompson et alj|2004) . but it poses a challenge to the hypothesis 
that magnetars are formed as a result of magnetic field amplifica- 
tion in a rapidly spinning (P ~ 1 ms) proto-neutron star. 

Note that when the magnetar theory of SGRs and AXPs was 
first proposed, the existence of magnetars was hardly an establ ished 
fact JDuncan & Thompsonll992tlThomDson & Duncanll993^ . To- 
day their existence is only doubted by a few, especially after the 
giant flare of SGR 1806-20. In general, SGR flare luminosities ex- 
ceed the Eddington luminosity, making it hard to find any other 
source of energy but the internal magnetic fi eld. The energ y of the 
giant flare of SGR 1806-20 was ~ 10'"' erg JHurlev et all2 005). If 
a neutron star has several of those flares this energy can only be pro- 
dded with an internal magnetic field of ^ 10^^ G (e.g. Stella et alJ 
l2005h . Moreover, the results presented do not rule out initial peri- 
ods of ~ 3 ms, the upper limit at which the a — Q. dynamo may 
operate JDuncan & Thompsonll 996). 

The generation of such a high internal magnetic field poses 
challenges in itself. The magnetar-sized magnetic fields may ei- 
ther derive from a high magnetic field in the core of the neutron 
star progenitor (the fossil field hypothesis), or may be generated in 



the proto-neutron star, in which case rapid rotation and strong con- 
vective motions and/or differential rotation are needed. It is indeed 
show n that proto -neutron stars are hig hly convec tive ( Buras et aj 
■2006:: iFrver &Wai-ren.2004.) . although lprver & Warrea 12004.) ex- 
press some doubt whether it is sufficient to generate magnetic fields 
in excess of 10^* G. The rapid rotation needed for amplifying mag- 
netic fields in proto-neutron stars may reflect the high angular mo- 
mentum of the progenitor's core. However, the proto-neutron star 
may also wind up due to asymmetric accretion during collapse, in 
which case it is expect ed that the accretion results in both addition al 
spin and kick velocity 'Spruit & Phiimevll998tlThomDsoij200Cl) . 

For the fossil field hypothesis, the solution to the question 
what the origin of the high magnetic field is, creates itself a new 
problem,namely, what is the origin of the high magnetic field of 
the progenitor star? It could be that the magnetic field is again the 
result of the amplification of a seed magnetic field due to rotation 
and convection, b ut in this case in the co r e of th e progenitor. How- 
ever, as noted bv lThompson & Duncanl il993ft . there is less con- 
vective energy available during any stage of stellar evolution, than 
during the proto-neutron star phase. Another option is what one 
might call the strong fossil field hypothesis; the magnetic field in 
the star reflects the magnetic field of the cloud from which the star 
was formed. 

Ironically, stellar evolution models indicate that progenitors 
with a high magnetic field in the core end up producing more slowly 
rotating pulsars I He ger et alj|2005l) . This is a result of the more 
effective coupling of angular momentum of the core to the star's 
outer radius, where angular momentum is lost due to stellar winds. 

The results presented h ere are in agreement with the fos- 
sil field hypothesis. Recently, iFerrario & Wickramasingha i2006h 
showed that the magnetic flux distribution of magnetic white 
dwarfs and neutron stars are similar Translated to neutron star 
radii, the highest magnetic field measured in white dwarfs corre- 
spond to 10"-10^^ G. However, this falls still short of the 10^^ G 
needed to power flares from SGRs such as SGR 1806-20. Possibly 
this could be accounted for by additional magnetic field amplifica- 
tion during the late stages of stellar evolution of stars that end in a 
core collapse (c.f. the core magnetic field evolution of the magnetic 
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star of lHeger et alllOOSi) . On the other hand, the energetics of the 
SNRs with magnetars could still be reconciled with magnetic field 
amplification in a rapidly rotating proto-neutron star, provided the 
angular momentum is not lost by magnetic braking, in which case 
it powers the supernova, but due to gravitational radiation or due to 
emission in a jet (see also Thompson & Duncan 1993) . Such a loss 
mechanism should work on a time scale shorter than the magnetic 
braking time scale, i.e. less than a few hundred seconds. 

Stella etal. (2005) argued recently that if internal magnetic 
fields exceed 5 x 10^*^ G neutron stars deformation is sufficiently 
strong to cause gravitational waves that radiate away most of the 
rotational energy. Note, however, that this is only true if rotation 
losses due gravitational waves exceed losses due to magnetic brak- 
ing. This requires Bi > 5 x 10^* G, for which no conclusive ev- 
idence exists yet, whereas the external dipole field has to be low, 
Bd < 5 X 10^* G. Moreover, this estimate is based on magnetic 
braking in vacuum, whereas during the first few seconds magnetic 
braking is likely to be enhanced due to the presence of hot plasma 
surrounding the proto-neutron star (Thompson et al. 2004). 

It is possible that in the early phase of a magnetar the dipole 
field is indeed low. On the other hand, the observed dipole fields are 
lower than is required to slow down to the spin periods of AXPs 
and SGRs. A case in point is lE1841-045/Kes73, for which we 
determined an age of ~ 1300 yr. To bring it from a period of few 
millisecond to its current period (11.8 s) a dipole field is required of 
1.6 X 10^^ G for a braking index of n = 2.5. Gravitational waves 
cannot explain this, since only for very short periods do gravita- 
tional wave losses dominates over magnetic braking, because ro- 
tation loss scales with Q:' for gravitational waves and with Q,"^ for 
magnetic braking. 

Another consequence is that as soon as magnetic braking dom- 
inates one expects the creation of a strong a pulsar wind nebula. In- 
triguingly, none of the magnetars studied here shows any evidence 
for a pulsar wind component. The radio map o f Kes 73 shows, i n 
fact, that the AXP is located in a hole in the SNR JKriss et alJl985^ . 
In X-rays we do not expect to see a pulsar wind nebula, as the 
ultra-relativistic electrons responsible for X-ray synchrotron emis- 
sion will have lost most of their energy. Note that the spin-down of 
a magnetar like 1E1841-045 is so rapid that the pulsar wind neb- 
ula must evolve during the early, free expansion phase of the SNR. 
Since most of the ejecta has not been shocked it is possible that the 
relativistic particles freely escape instead of being confined by the 
hot plasma of the SNR. However, it would be of great importance 
to investigate the formation of pulsar wind nebulae from magnetars 
in more detail, as it may provide additional means of determining 
whether AXP/SGRs started with periods in the milli-second range. 

Finally, another way of losing angular momentum is by means 
of a jet during the explosion. Radio and X-ray maps of Kes 73 
and N49 do not show evidence that a jet may have been present in 
the past. Kes 73 is in fact remarkably circular symmetric (Fig.0. 
Interestingly, however, there is considerable evidence for a jet in 
the SNR Gas A (Vink 2004; H wang et al. 2004) . and there is indeed 
evidence, al though no t conclusive, that the point source in Gas A is 
a magnetar JKrause et ali2005i) . 
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Table Al. Best fit spectral models to XMM-Ne^vtonlMOSl data. 





Kes73 


N49 


Model 


2NE1 


Sedov 


Sedov 


EM(10l2cm-5) 


26.9 ±3.1 


16.4 ±1.5 


5.4 ±0.5 


fcTe (keV) 


0.63 ±0.06 


0.72 ± 0.3 


0.415 ±0.005 


net (10"cm-3s) 


3.1 ±0.8 


4.6 ±0.2 


28.0 ±2.6 


EM2 (10i2cin-5) 


2.0 ±0.5 






fcTao (keV) 


2.26 ±0.35 






net2 (lO^cm-^^s) 


0.47 ±0.10 









(1) 


(1) 


0.25 ±0.02 


Ne 


(1) 


(1) 


0.41 ± 0.02 


Mg 


1.13 ±0.13 


1.95 ±0.15 


0.27 ±0.01 


Si 


1.09 ±0.07 


1.7±0.1 


0.45 ±0.01 


S 


1.12 ±0.06 


1.6 ±0.1 


0.99 ±0.06 


At 


1.02 ±0.18 


1.1±0.1 


- 


Ca 


1.84 ±0.42 


0.77 ±0.33 


- 


Fe 


0.42 ±0.15 


4.3 ±0.7 


0.19 ±0.01 


ArH(102icm-2) 


27.3 ±0.6 


31.2 ±0.6 


2.2 ±0.1 


Fit range (keV) 


0.8-8.0 


0.8-8.0 


0.5-3.0 


C-statistic/d.o.f. 


186/110 


956/468 


574/153 



NOTE - Errors correspond to 68% confidence ranges (AC = 1.0). EM 
refers to the Emission Measure ( / ricnudV / {JTzd?)). Abundances are 
given with respect to the solar abundances of lAnders & Grevessd 11989B . 
For Kes 73 the abundance of oxygen and neon were fixed to solar values, 
because the strong absorption did not allow for an accurate determination. 

APPENDIX A: THE SPECTRAL ANALYSIS OF 
XMM-NEWTON DATA OF KES 73 AND N49 

For determining the energies of Kes 73 and N49 we used archival 
XMM-Newton data observations made on October 5 and 7 2002 
(Kes 73,ObsIDs 0013340101 & 0013340201) and of N49 made 
on July 7 2000 and April 8 2001 (N49, ObsIDs 0111130301 & 
01 13000201). We only used data from the MOS 1 and MOS2 instru- 
ments , because they offer the best energy resolution JTumer et alj 
I2OOII) . We processed the data in a standard manner employing the 
XMM-Newton software SAS 6.5.0, excluding time intervals with 
excessive background flares. The source spectra were extracted 
from circular regions encompassing the source, excluding small 
circular regions containing the AXP or SGR. Background spectra 
were extracted from annular regions around the SNRs. For our fi- 
nal spectral analysis we added all the MOSl and M0S2 spectra 
of each individual SNR together for statistical reasons, using expo- 
sure weighted merged response files (RMFs) and ancillary response 
files (ARFs), a procedure that makes for a nicer representation of 
the spectra and gives little loss of accuracy, as the two MOS detec- 
tors are similar. The total effective exposure per detector is 9.8 ks 
for Kes 73 and 39.6 ks for N49. 

For the spectr al analysis two spectral analysis packa ges were 
employed, SPEX ^ Kaastra & Mewdl2000 h and XSPEC JAniau^ 
[1996). SPEX has the advantage that it includes more lines as- 
sociated with inner shell ionization, which is important for non- 
equilibration ionization (NET) plasmas. This is in particular rel- 
evant for Kes 73, which is more out of ionization equilibrium 
than N49^ (Table |aIJ. XSPEC contains a Sedov model (vsedov, 
[Borkowski et al. 200ii) that also includes non-equilibration of elec- 
tron and ion temperatures, an effect that is known to be important 



for in particular SNRs with a s hort ionization age (rict) and high 
shock velocity ( Rakowski 2005). We used SPEX for a more heuris- 
tic approach in which we fit the spectra with two non-equilibration 
ionization models with two different temperatures and ionization 
ages, but with elemental abundances of the two components cou- 
pled. 

From the results listed in Table IAII one can calculate the pre- 
shock density using the relation for SNRs in the Sedov phase 
J rianudV — 0.129^rJ. The shock velocity can be obtained 
from Eq. |3| and using the relation between Vs , r^ and t (Eq. 
one can then infer the age of the SNR. From the age and the 
pre-shock density an estimate of the explosion energy can be de- 
rived (Eq. 0. However, Eq. |3|does not apply in case the plasma 
is not in full temperature equilibration. Although both remnants 
have lower shock velocities and much larger ionization ages than a 
SNR like SN1006 for which this effect h as in fact been measured 
(nct= 2 X 10'' cm'^s. l^iik et alJl2003h . there may still be some 
effect, in particular for Kes 73. The plasma parameters of Kes 73 
are in fact quite similar to those of CTB 109, for which Sasaki et alj 
1 2004) argued that Tp/Te ^ 2.5. We find that for Kes 73 assuming 
non-equilibration gives a more consistent result. The reason is that 
the age of the remnant can be estimated both from kT and from 
the ionization age. From rict we derive t — 1300 ± 200 yr, but 
using kT we find from the values in the second column of Table Q 
2145±37 yr for full equilibration, but 1357±23 yr if Tp/Te = 2.5. 
For N49 the ionization age and temperature give consistent results 
assuming full equilibration (resp. 6.6 ± 0.9 kyr and 6.1 ± 1.7 kyr). 
This is to be expected, since the ionization age of N49 corresponds 
roughly to the Coulomb equilibration time. 

A note of caution concerning the interpretation of the X-ray 
spectrum of Kes 73: As Table lAll shows. there is considerable un- 
certainty about the abundance pattern of Kes 73, with the 2NEI 
model suggesting abundances that are almost solar, whereas the Se- 
dov model suggest an enriched plasma. Other differences between 
the models in fact go back to this, because fitting the Sedov model 
with the abundances fixed at solar values we obtained a total emis- 
sion measure and temperature closer to that of the 2NEI model. 
The matter is of some interest, as overabundances would be an in- 
dication that Kes 73 is in an early evolutionary state that cannot be 
well described by a Sedov model. Another possible indication for 
that is the low total mass of Kes 73 (25-33 Mq, for resp. the 2NEI 
and Sedov model), much lo wer th an N49 (320 Mq) or CTB 109 
(97 ± 23 M,7.. ISasakietalJl2004h . We consider the 2 NEI model 
to be more reliable in this case, since SPEX incorporates more line 
data relevant for NEI plasmas. In fact. Fig. Q shows that the Fe-K 
complex around 6.6 keV and the Mg XII lines around 1.4 keV are 
not well fitted. The SPEX 2NEI model does not have these prob- 
lems. Although, the fits to the spectra of Kes 73 and N49 are not 
perfect, the global features of the spectra are well fitted (Fig.0. 

Nevertheless, we tried to constrain the explosion energy of 
Kes 73 further by exploring the more generic SNR evolution mod- 
els of Truelove & McKee (JW99). We tried to answer the ques- 
tion, which explosion energies are still consistent with the densi- 
ties and post-shock temperatures found from our fit, 2-4 cm~"^ and 
0.6—1.9 keV, given the shock radius of 4 pc and assuming that the 
ejecta masses lie in the range 0.5-25 Mq. Considering ejecta den- 
sity profiles with power i ndices of n = 0, n = 7, and n — 9 (see 
iTruelove & McKeell999l) . it turns out that the observational data is 
only consistent with energies in the range (0.1 — 1.1) x 10^^ erg. 



t For a recent detailed X-ray study of this SNR see lPark et alj l2003h 
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